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1 .0 INTRODUCTION 


Ihis third monthly progress report contains corrections and additions 
to previously submitted reports. The additions include the results of 
regression analyses on the tabular aerodynamic data provided. The objectives 
of the regression analysis were two-fold. The first was the establishment 
of a representative aerodynamic model for future coefficient estimation. 

The second was to reduce the storage requirements for the "nominal" model 
used to check out the estimation algorithms. The results of the regression 
analyses are presented in l^pendix D. 

The primary activities during this reporting period were the develop- 
ment of the computer routines for the filter portion of the estimation 
algorithm and the "bringing-up" of the NASA SRB predictive prc^gram on SDI's 
computer. For the filter program i approximately 54 routines have been 
developed. The routines have been highly subsegmented to facilitate over- 
laying program segments within the partitioned storage space on our computer. 
The attempt to bring up the NASA SRB predictive program on our computer has 
not been completely successful to date. With such a large program, over- 
laying techniques are required. Progress has been made and we are 
"hopefully" in the final stages of shifting subsegments around within the 
segments such that no segment is larger than the allocated partition. 


2.0 FILTERING AND SMOOTHING ALGORITHM 


The Space Shuttle Parameter Estimation Program utili.zes optimal 
estimation techniques to provide estimates of the propulsion system 
parameters . The technique selected is the extended Kalman filter and 
the modified Bryson -Frazier smoother. By modeling the propulsion system 
parameters as time correlated random variables » improved estimates of 
these parameters are obtained and are properly time phased by removing 
the filter induced lag by using the combined fil ter /smoother . The 
smoother also provides impr 'ved estimates of the initial state estimates 
The system, in state-space notation, is modeled as the continuous 
dynamical system disturbed by additive Gaussian white noise 

X = ^(x(t) , t) + G(t) w(t) + u (t) , x(o) ~ ^ (1 ) 

where 

X = n -dimensional state vector 

X = Gaussian initial condition vector with covariance P 
—o o 

w(t) = p-dimensional vdiite, zero-mean white Gaussian noise with 
covariance 

E[w(t) w’^(t)] = Q(t) 6(t - t) 
u{t) = n-dimensional control vector. 

The elements of the vector x(t) represent vehicle position, velocities 
attitudes , angular rates , aerodynamic and propulsion parameters , measure 
ment biases, etCs Elements of u(t) include known control inputs such 
as SSME power level commands . 
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The system described by equation (1) is observed at discrete times > t 
with not all states being directly measured. Some measurements are non- 
linear functions of the elements of the state vector general 

the measurement process is described as 

\ ^2) 

where 


z, » m-dimensional observation 
-k 

K - functional representation 
"nC 

the states 

V, = m-dimensional > zero-mean • 
covariance 


vector 

of the measurements in terms of 
while Gaussian noise sequence with 


E[v. vT] t= R. 
“1 -“3 1 


6 . . 
3 


Examples of the elements of the observation vector ^ include radar 
measurements of range, azimuth, and elevation from the radar site to the 
vehicle . 

It is assumed that the system process noise vector w(t) and the 
measurement noise vector are uncorrelated. Also, the system state 

initial condition vector x is not correlated with either of these two 

mmQ 

noise vectors. Therefore 

E[w(t) = 0, E[w(t) x^] = 0, E[x = 0 

~K "*0 "TC 


where the superscript T denotes transpose. For later reference, the 
following matricies are defined 


OF pS ^ 
POOR QUALiry 


F(x(t), t) 


9f^(X(t), t) 
0x(t) 


(3) 


and 


H(5^{-)) 


8h(x(tj^) ) 
3x(t^) 


(4) 


These matricies are linearizations of the dynamics and measurement models 
respectively i evaluated about either a nominal or reference value of the 
state) or about the state estimate. 


2 . 1 Extended Kalman Filter Algorithm 

The extended Kalman filter algorithm is in essence a conventional 
linear Kalman filter algorithm applied to a mathematical model resulting 
from the linearization of the system model equation (1)) and measurement 
process) equation (2), about a current state estimate w The filter yields 
optimal estimates if the linearization is accurate) i„e.) the state esti- 
mate closely approximates the true state. The derivation of the algorithm 
can be found in reference [1 ] . 

The algorithm proceeds as follows. After initialization of the state 
estimate and covariance) tha state estimate and covariance are propagated 
forward in time until a measurement update is available) by 

X = f.(x(t)) t) , tj^ ^ < t £ tj^ (5) 

and 

P(t) = F(x(t)) t) P(t) + P(t) F(x(t)) t)"^ + G(t) Q(t) G(t)^ 


( 6 ) 
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At the measurement timoi the state estimate and covariance are updated 
by 


^{ + ) = ^(-) + - h^(^(-))) (7) 

and 

P^( + ) = (I " \ 

where the (-) and (+) represent the appropriate values just before 
and just after the update. The updated values are used to reinitialize 
the time propagation equations (3) and (4) for integrating up to the next 
measurement time . The Kalman gain matrix is computed as 

This algorithm is repeated until the last time point, t^^, is processed. 

For later use in the smoother algorithm, various combinations of the state 
estimates l,x), measurements (^), linearized dynamics matrix (F) and 
measurement matrix (H), measurement noise covariance (R) and estimation 
error covariance matrix (P) must be stored for each time instant to be 
processed by the smoother algorithm. 

2 .2 Modified Bryson^-Frazier anoother Algorithm 

The operation of the smoother algorittim is similar to the filter 
algorithm except in reverse time. The derivation of this smoother algorithm 
is found in reference [2]. This fixed interval smoothing algorithm pro“ 
vides optimal estimates given all the measurements in comparison to the 
filtering algorithm providing optimal estimates given the previous 
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measuremonts procossod. Therefore the smoother provides improved estimates 
in addition to removing the time lag induced by the filter algorithm. 

The smoothing algorithm adjoint variables , ^ and A are "initialized" 
at the final time processed by the filter, T, 

XW-) = -h;;(h„ h;; + . h„(s^<-))) 6 do 

N , T 

and 




N,T 


(11 ) 


If T is not an observation time, ^ and A are zero. The adjoint 
variables are propagated in reverse time to the next previous measurement 
time by 


X = - P(x(t), t) X , ^ ^ \+1 

A = - F(x(t), A - A F(x(t), t) (13) 


At the time of an available measurement, t, , the adjoint variables 


are updated by 

X(-) = X(+) - hJ(H^ ^k ^ 

and 

A(-) = (I - Kj^ A( + ) (I - 


( 14 ) 
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The smoother state estimate and error covariance are obtained using the 
filter estimate and covariance and the adjoint variables by 

x*(t) « x(t) ~ P(t) X(t) (16) 

and 

P*(t) = P(t) - P(t) A(t) P(t). (17) 

Due to the potential number of time points to be processed, smoother 
estimates may only be computed at the discrete measurement times. For 
this approach the propagation equations (10) and (11) are replaced by 

and 

T 

where rs the state transition matrix formed with the linearized 

dynamics matrix F to propagate the adjoint v'ariable from time 
to time tj^. The algorithm continues in reverse time until the initial 
time is reached. 

2 .3 Iterations with the Filter/anoother Algorithm 

The performance of the filter /smoother algorithm is a direct result 
of the accuracy of the linearization. Repeated operations of the algorithms 
with adjustments in initial state estimates and covariance in each cycle 
can yield improved estimates . This technique is known as global iterated 


I 



tiltoring as dofinod in roforonco [3] . Each cycle of operating the 
algorithms would yield increasing improvomentB in the state estimates. 

This feature of the algorithm fffXiration is of special interest to 
the propulsion parameter estimation problem using the NASA predictive 
models. Initial, or nominal, values of the parameters of interest can 
bo used to obtain the necessary partial derivatives indicated earlier. 
Prom operating the algorithm improved estimates of those parameters are 
obtained. Using these ojnproved estimates, more accurate partial deriva- 
tives are obtained for us© in the algorithms. This process is continued 
until there is in essence no change in the partial derivatives or quality 
of the state estimates. If the linearization is accurate, the measure- 
ment residual should be a white noise process with known covariance . 


3 . 0 FILTER/SMOOTHER ALGORITHM SYSTEM AND MBASUREMEN'f MODEL 


Iho UBofulnosc of the filter /smoother algorithm is to provide esti- 
mates of the system states from the observed motion and dynamics while 
the systt^m is driven by known and unknown elements. These unknown elements 
are elements of the system state vector to be estimated. The evolution of 
motion resulting from those known and unknown elements is assumed to bo 
suitably represented for this study by a six degree-of -freedom (6 DOF) rigid 
body equations of motion. These equations are presented and discussed in 
section 3.1. 

To implement these equations into the filter /smoother algorithm 
presented in section 2.0, a linearization of the system state and measure- 
ment models is required. These linearized equations are presented in 
section 3.2. 


3 .1 Equations of Motion and Measurement Equations 


3.1.1 Rigid Body Equations of Motion 

( B ) 

The rate of change of vehicle velocity in body coordinates , v > 


as a result of external forces acting on the vehicle is described by 
2 


(B) P^'^m 


2m ^ 


c. . (r<^^ 


f(B) .(B) 

(B) -T ~P 
U) X v + - 


m 


m 


( 20 ) 


where 

p = atmospheric density 

A = aerodynamic coefficient referenced area 

v^ = magnitude of vehicle velocity relative to the surrounding air 


mass 
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m c vohxclo mass 

B aerodynamic force coefficient vector 
(I ) (I ) 

2. (r ) » gravity vector in inorca.al coordinates 
^ B angular rotation of the body relative to the inertial frame 

B resultant thrust force vector in body coordinates 

f « resultant plxme force vector in body coordinates 

( 1 ) 

llie rate of change of vehicle position in inertial coordinates, £ > 

is then obtained by 


(!) 




(21) 


X B 

where C is the transformation matrix from body coordinates to 

T Q 

inertial coordinates . The elements of the C transformation matrix 
are obtained from the resulting Euler angles defined by 



• 

<P 


1 

sin<ptan0 

coseptanO 


P 


• 

e 

= 

0 

COSfp 

-sin<p 


q 




0 

sin<psec0 

costpsecO 


r 

where 

9 


, and 

<j> are roll. 

pitch and 

yaw 


The roll, pitch and yaw rates of the body relative to inertial coordinates 
are p, q, and r respectively. Finally, the rate of change of the body 
rates relative to inertial is given by 
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UJ 


IP 


i q 


ur' [ 


PAv" 


m 


Phv 


“-~tn 


m 




-eg 


X 


- W X (TO)) 


ir 




(23) 


whoro 


I b: vehicle moments of inertia matrix 

Zc » aerodynamic coefficient referenced length and moment 
'""■’fn 

coefficient vector 

r^®^ » vehicle centcif'*of -gravity vector in body coordinates 

—erf 

■*"13 

* t*' ^ 

i ' B aerodynamic coefficient reference position in body coordinates 
-“A 

/ n \ 

^ » resultant thrust torque vector in body coordinates 
b resultant plane torque vector in body coordinates 


The equations of motion represent the first twelve elements of the system 
state vector. These equations are summarized in Table 3.1 ,1-1. 

The moment of inertia matrix I in general is given by 



-I 

-I 

X 

xy 

zx 


I 

-I 

xy 

y 

yz 


-I 

I 

zx 

yz 

z 


L J 


(24) 


for the moment axis terms i i.e., , and the product of inertia terms i 

i . e . , I 

zx 
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The aorodynamic t'orco and moment coefficients and plume forces are 
defined as functions of angle -of -attack, ot, and angle-of -sideslip, 3, as 
shown in Figure 3.1 .1-1. The body referenced relative velocity vector, 
removing the wind velocity, v , from the vehicle velocity, is given by 


(B) B^I (B) B_LL (LL) 

Vi=V -CV»=V - C V 
~T — — W — 


(25) 


where v 


(LL) 


is the local -level referenced wind velocity vector. Ihe 


following equations define a and 3 in terms of the components of v 


V 

a ss tan"^ 


(26) 


sin (- 


m 


where 


(27) 


, 2,2 2 , 

V =: (V + V + V ) 

m r^ r^ r 3 


(28) 


(b ^ 

The resultant thrust force ' is expanded as 


“T 1=1 1 


f^ 

T. 

1 

0 

0 


“ iii *^i 


(29) 


B ^ 

where the transformation matrix transforms the magnitude of thrust 

for each thrusting device, f,^, , from its center -line to the body 

■^i 

coordinates . The general eqaution for f^ is 

i 



UJ>. 

« t ce 

HOO 

<o»- 

J U U 
UJ lU tu 

ts>> 



FIGURE 3 .1 .1 -1 . Body Referenced Wind Axis System 



-PA 

T. T. s. e 

X xvac X 


where 


= vacuum thrust 

ivac 


P = atmospheric pressure at motor exit 
“i 


A = motor exit cone area 
e 


B ^ 

The matrix is different for the SSME's and SRB’s and is given by 


n 

B X 
C 


B MP MP G, G 3 


B ^ 
C 


SSME 


SRB 


( 30 ) 


where 


B MP 

C = transformation from the engine mount plane to the body 
coordinates 


MP_G 


C = transformation from the gimbal reference plane to mount plane 
(structural deformation) 




C = transformation from enterline to the gimbal reference plane 




C = transformation from SRB nozzle centerline to the body 
coordinates (gimbal angles ) . 
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The resultant thrust torque is the summation of the torque contribution 
from each thrusting device and is given by 




where 


iSi 


'ijp 

i 


(B). 

r ) X C. 
-eg 1 


■T. 

1 

0 

0 


(31 ) 


r(B) 

1 


body coordinates of the thrust reference point for the 
i^^ thrusting device. 


3.1 .2 Measurement Equations 

The measurements assumed available for the filter /smoother algorithm 
include inertial platform acceleration and attitudes j ground based radar 
tracking, SRB' s head pressure, SSME' s chamber pressures, liquid flow 
rates, pressurant flow rates. The ET volumetric levels are available; 
hov/ever, due to their limited nmber (4), they may only be used for 
alternate checks of the filter /smoother algorithm performance . 

The propulsion related measurements will be treated in a separate 
section. In the following, the inertial platform acceleration measurements, 
attitude measurements and ground based tracking measurements models will 
be described for later linearization. 
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3.1 .2.1 Platform Accoleratlon Measurements 

Accelerometers mounted orthogonally on an inertially stabilized 
platform, not located at the vehicle center of gravity, sense externally 
applied special forces and accelerations due to body rotation, Ihe 
accelerometer measurement is modeled by 




S^P P^P* 

= c c 


P-CB 


pAy 


m 


2m 




.(B) (B) 

Lp 

+ 


m 


m 


+ W X (1) X 




-eg 


+ W X 




-si!’” 

-eg 


4“ D 4* V 

—a 


(32 


where 


S P 
C 

P P' 

c 


P' B 
C 



b<®> 

—a 


V 
— Q, 


(S) 


= transformation from platform coordinates to sensing coordinates 

= transformation from misaligned platform coordinates to 
platform coordinates 

= transformation from body to misaligned platform coordinates 
= body coordinates of the platform center 
= accelerometer bias vector 
= accelerometer measurement noise vector 


3.1 .2.2 Platform Attitude Measurements 


The inertially stabilized platform for the STS is a four axis IMU 
with a redundant roll axis [4] . Vehicle body attitiides are generated 
via quaternions [5] . It is assumed that an equivalent representation 
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can be made to obtain vehicle attitude by a three rotation sequence of 
roll, pitch, yaw to transform from inertial to body coordi antes. Hiis 
approach has been used in reference [6] . 

The attitude angle measurement model is given by 


6 <=> 

—m 


0 + 




(S) 




(S) 


where 


(33) 


P P' 

^ = platform misalignment bias vector (used to formulate C ) 




( s ) 

' = attitude measurement noise vector. 


The transformation matrix used to trcinsform from body to inertial 
coordinates in terms of the elements of the 6 vector is given by 


cos6cos{|j sin(psineoos(jj 

-cos(psin(l> 


cos<psin9cos(j) 

+sincpsin(|j 




cos0sin<jj 


sin(psin0sin(p cosq>sin0sin(l< 

+Gos(pcos(jj -sin<pcos«}j 


-sin0 


sincpcosO 


costpcos0 


(34) 


3. 1.2. 3 Ground Based Tracking Measurements 

Ground based radar tracking devices can provide measurements of 
range, azimuth and elevation from the radar sight to the vehicle. Azimuth 


and elevation are established relative to the sight's local level. If 
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tho txacking device is a passive optical tracker (not laser) then only 
azimuth and elevation measurements are available requiring more than one 
to establish position information. 

Defining x, y, and z as the local easti north and up position of the 
vehicle relative to tho ground based tracking device, the radar measurement 
equations are given by 


/ 2 , 2 2 / . 

p=(x +y +z) +b 

P 

A — tan ^ (~) + b + v 

V A ^ A 

E = tan"^ (- - ^ ) 

r~2 ~ 

y X + y 


+ V 

P 


+ b_ + AE + v_ 

tj hi 


( 35 ) 

( 36 ) 

( 37 ) 


where 


b , b , b_ = range, azimuth, elevation biases 
p A E 

AE = atmospheric refraction correction 

Vp, v^, Vg = range, azimuth, elevation measurement noise. 

•The position vector of the vehicle relative to the tracking device is 
given by 



y 


A . (LL) LL^ECP rECP^ECI (I) (ECF), 
= toy = c [ c r - rgpR ) 


( 38 ) 


z 


whore 
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LtL ECP 

C e transformation from earth center fixed to local level 
ECP ECX 

■' C = earth centered inertial to earth centered fixed 
/ ECP ) 

“ positlon voctor of tracking device in ECF coordinates . 

^\UK 

LiLi ECP 

The transformation matrix C is given by 


sinX 

-sinLcosX 

cosLcosX 

cosX 

-sinLsinX 

cosLsinX 

0 

cosL 

sinb 


(39) 


where L and X are the geodetic latitude and east longitude of the 

ECP ECX 

device. The transformation matrix C is given by 


ECF^ECI 


cos [to (t — t ) ] 

-sin[o) (t - t )] 
e' RNP 


oos[u>p,(t - tp^p)) 


[RNP] (40) 


where 


0 ) 5= earth rotation rate 

E 

t_„_ = time tag for RNP matrix 
RNP 


The position vector , r„__ , 

“*"KDR 


of the tracking device is given by 
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(ECP) 

-RDR 


a 


— + h) eosLcosX 

V cos^L + (1 ~ e)^ sin^L 

( — - - + h) cosLsinX 

\l 2- v2 . 2 

V cos L + (1 - e) sin L 

R„(1 - e)^ 

( — - ’ h) sinL 

\f 2, .2 . 2, 

\ cos L + (1 - e) Bin L 


whore 


R„ =s equatorial radius of Fisher ellipsoid 

Cj 

e a flattening of Fisher ellipsoid 
h = altitude of the device above Fisher ellipsoid 


3.2 Linearized System State and Moaouromont Equations 

The vehicle equations of motion are nonlinear functions of their 
motion variables and are implicit functions of other elements of the 
system states . The measurement equations involve similar function rela- 
tionships. Ihe linearizations for the filter /smoother algoritlim require 
partial derivatives with respect to the motion variables, i.e., v ^ and 
0, and with respect to other elements of the state vector, yielding 
explicit functional relationships for the elements of interest. 

For system state equations the partial derivatives will be presented 
in section 3.2.1 for the state elements in order of occurrence for the 
first twelve states. Other partial derivatives for candidate state ele- 
ments will follow in section 3.3.1. ihe measurement equation partial 
derivatives for the first twelve states will be presented in section 3.2.2. 
Partial derivatives of the measurement equations for other candidate states 
will bo presented in section 3 .3 .2 . 

Hie resulting partial derivatives are imbedded into the linearized 
system state matrix, P(x(t), t), as shown in Figure 3.2-1 . A corresponding 

linearized measurement matrix, H(^)^ is similarly formed with the 
measurement equations' partial derivatives. 

3 . 2*. 1 System State Partial Derivatives 

Peirtial derivatives of each of the equations listed in Figure 3.2-1 
are developed in their order of occurrence with respect to the order of 
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Inertial Poaition Rato Equation 

The first nonzero partial derivative of tho r^^^ equation is with 
respect to 


av 


(B) 




I B 

C . 


{42) 


Tho second nonzero partial derivative is with respect to 0 . This partial 
derivative results in a third order tensor and occurs frequently in later 
developments. The generalized form is presented in Appendix A. 


Body Velocity Rate Equation 

• (b ) ( I ) 

The partial derivative of v' with respect to r is given for 


altitude terms approximately as 

,j,T 

>«* 

a 


dr 


(I) 


d_ ~ 

ah , (I) 


(43) 


where 
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1110 partial dorivativos of 
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and 


m 
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occur frequently and are given in Appendix B. 

B I (I ) (I ) 

The gravity vector C (r ) partial derivative with roopect 
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where 


U *3 gravitational constant. 

The partial derivative , — rry > is the suin of the matricies in equations 

Or' ' 

43 and 45. 

p f R ^ ( R 5 

The partial derivative of v' with respect to v' ' is given by 
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whore 


B skew symmetric matrix made from the elements of the vector m 
and equivalent to the cross product operator w x ( ) . 

• (B) 

Oiie partial derivative of v wxth respect to ^ is 
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The partial derivatives of are given in Appendix B and the partial 


-r 


3v 
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derivative — is given in Appendix A. The last partial derivative 


is given in Appendix C. 
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The partial derivative of v with respect to w is 
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Eulor Angle Rate Equation 


The Euler angle rate equation is a function of both the Euler angles 
and the inertial rates . The linearization will yield the two associated 
matricies . 

First with respect to the vector 6_, the following matrix results 


ae 
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The partial derivative of 


0 with respect to w is 
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Inertial Angular Acceleration Equation 

The first partial derivative of this equation is with respect to the 

(I ) . 

vector r . Using the approximation indicated in equation 43 , this 


partial derivative is 
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Hie partial derivative with respect to the vector 0^ ia 
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The final partial derivative for the first twelve states is with respect 
to the vector u. This operation yields 
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3.2.2 Measurement Partial Derivatives 

The measurements assumed to be available, as discussed earlier, include 
ground based radar tracking, inertially stabilized platform attitudes 
relative to the vehicle body, and stabilized platform mounted 3 axis 
orthogonal accelerations. As with the state dynamics matrix, the measure- 
ment equations are linearized about the best state estimates. 


Radar Track Measurement Equation 

Referring to the radar track measurement equations, the required 
partial derivatives are 


i!£. 


8r 




—V 


(I) 


-V 


8r 


(I) 


8r 


8A 

(I) 


8A 


aAr<^^> 


-V 


8r 


(I) 


8r 


8E 

(I) 


8E 


6Ar<“> 

—v 


aAr<^^' 

—V 


6r 


(I) 


(55) 


(56) 


(57) 


The last partial derivative’ in each of these equations, 


aAr<^^> 

—V 



is 


8Ar 

—V 


8r 


(I) 


LL ECF ECP^ECI 
C C 


(58) 



ORfGJWAL PfC2 
OF POOR quai;t\" 


The rest of the recjuired partial derivatives are 
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Inertially Stabilized Platform Attitude Equation 

The inertial platform is assumed to provide attitude angle measure- 
ments of the true attitude plus an attitude bias plus measurment noise . 

The partial derivative of the measured attitudes with respect to the 
vector 0 yields an identity matrix. 

Accelerometer Measurement Equation 

The accelerometer senses specific body forces excluding gravity along 
the sensing axes. With reference to the accelerometer equation^ the 

(I) 

partial derivative with respect to r is 
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The partial derivative with respect to v' yields 
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For the partial derivative of the accelerometer with respect to 
the vector 0, the measurement equation is temporarily rewritten as 


^(S) S P' P' B _^(B) , , (S) ^ , (S) 

a = C Cs 4-t) +v 
— m — —a ~a 


where the vector s_^ represents the sum of the aerodynamic, thrust, 

pi B 

plume and rotational coupling terms. Ihe matrix C is the same 


matrix as C . The required partial derivative results from 
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“al7 = ^ 90 C S . 


( 65 ) 


The partial derivative on the right hand side is developed in Appendix A 

/ n \ 

with the vector s' ' representing the sum of the terms indicated above. 

The final partial derivative for the accelerometer measurement is 
with respect to the body rotation vector m. Defining 
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3.2.3 Additional Parameter Partial Derivatives 

The mathematical developments are presented in this section for the 

partial derivatives of the system and measurement equations to allow for 

additional candidate parameters to be included in the estimation algorithms. 

Those parameters include center-of«-gravity, £^g» moments of inertia, I, 

wind velocity, v , and inertial platform tilt errors. Aerodynamic and 
*"*W 

plume parameter partial derivatives are also presented. 

The computer program is being structured to permit these parameters 
to be easily incorporated without significant impact on the program code. 


3. 2. 3.1 Center-of -Gravity 

From equation 23, the partial derivative of angular acceleration with 

respect to r is 
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From equation 32 , the partial derivative of the measured acceleration with 

respect to r is 
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3. 2. 3, 2 Moments of Inertia 


The moments-of -inertia are grouped into "principal" terms, i , and 

“P 

cross product terms, i . From equation 24, these vectors are defined as 
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With these definitions, equation 23 is rewritten as 
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where represents the sum of the nonrototional torques in equation 23. 

Defining an intermediate vector a as 

a_ c ST *■ w X (iw) (73 ) 

the partial derivatives of the angular acceleration with respect to 

i and i are 
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3. 2. 3. 3 Wind Velocity 

From equations 20 and 25 » the partial derivative of the vehicle 
acceleration with respect to v is 
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The first of the partial derivatives in equation 79 can be obtained 
from the following equation 
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From equation 80 > the following is obtained 
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Denoting tho olomonts of tho matrix \ as the following 

equations are obtained for tho partial dorivativos of « ond & with 
rospoct to V ; 
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3. 2, 3. 4 Inertial Platform Tilt 

Temporarily rewriting equation (32) as 
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6_Q = vector whose elements are the axes misalignments 

P' B 

s «3 sum of the bracketed terms in equation 32 multiplied by C . 


The following partial derivative of the measured acceleration with 
respect to 60 is obtained 
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3. 2. 3. 5 Aerodynamic and Plume Parameters 

A linear model for the aerodynamic and plume characteristics is 
used. This model is expanded as 

£f = ® 


and 


c = c + c a 
-m ~mo —mot 




(87) 


f=f -t-f a+f-3 + ... 

-p -p -pa ~pp 


( 88 ) 


where additional terms to represent rates, cross couplings, and controls 
can be included. 

The basic approach of establishing the partial derivatives will be 
illustrated for a couple of terms, c^^ and • Using these example 
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illustrations > the rest of tho candidate parameters can be similarly 
obtained. Prom equation 20, the following partial derivative is obtained 
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where 


[U] = unit 3x3 matrix with one's (1) on the diagonal and zeros 
off the diagonal 


Prom equation 23, the partial derivative ef angular acceleration with 

respect to c is 
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The corresponding partial derivative with respect to c^ is 
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The static aerodynamic coefficient model has been obtained by a multiple 
regression analysis of the current aerodynamic tubular data. This model 
is presented in Appendix D with the associated regression coefficients. 


3 .3 Propulsion Paramotor States and Moasurements 


A candidate approach for incorporating the NASA propulsion model's 
capabilities has been identified. This approach utilizes nominal pre- 
dicted values of thrust, pressure, propellant and pressurant mass flow 
rates, and utilizes sensitivities or partial derivatives of these 
variables with respect to the independent parameters selected for esti- 
mation by the algorithm . 

The approach is to include deviations from nominal values of 
measured chamber pressure, povrer level, propellant and pressurant mass 
flow rates as states. The models assumed for these deviations are time 
correlated random processes. Then as states, partial derivatives of 
the first twelve states with respect to these variables will be required. 

For the SSME and SRB, this modeling approach is discussed in the 
following. Mditionally, the necessary partial derivatives of the first 
twelve state variables with resepct to the additional states are presented. 

3 .3 .1 SSME Propulsion Parameter Model 

For the SSME, the total actual values of vacuum thrust and oxidizer 
mass flow rates are modeled by 

^t“=^t 

nom 

and 

mg = m^ + Airig . (93) 


IhQ measuremGnts of fuel mass flow rate, pronsurant mass flow rates and 
power level are modeled as 
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Ihese measured quantities include measurement noise s ^ ^ and potential 
bias states b^ ^ modeled as random constants. In these measurements > 
the A'd variables are to be included as states in the estimation 
algorithm. If the nominal values are zero or unknown, then the A'd 
variables absorb the entire estimate. Where required, the estimate for 
the variables used in the estimation algorithm is foirmed using the 
nominal and the estimate of the deviation, etc. In example, thrust and 
fuel mass flow rate estimates are formed as 
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and 
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The deviation or A'd measurement variables are modeled as time 
correlated random variables . This permits these veuriables to vary within 
a band of frequencies. The typical model is then given as 

A( ) = - ~ A( ) + 1 B. . (100) 

where the parenthesis ( ) would be replaced by the variables, i.e., m„ . 

* 2 
For the SSME, an additional variable modeled as in equation 

100 and included as a state variable with the A'd measurement 

variables . 

The thrust deviation is expanded as in the following truncated Taylor 
series as a function of the independent parameters. 
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In the V 
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and (I) equations, with equation 101 replacing f 
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the partial derivatives of f^ with respect to the A'd variables are 
obtained directly from equation 101. 
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It is desirable to include vehicle mass bias as a state. Ihe SSME's 
system contribution to the mass deviation is given by 
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In equation 101, the Am^ contribution to the mass deviation is not 
available from measurements. As with the thrust deviation j this quantity 
is formed as 
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which is in terms of other estimated state variables. In equations 101 
and 103 the deviation in mixture ratio, AMR, is obtained algebraically 
from 
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3.3.2 SRB Propulsion Parameter Model 


The approach for the SRB modeling follows closely that used for the 
SSME. Candidate independent parameters include propellant burn rate 
exponentj a, and motor efficiency coefficienti . Others can be added 
using this technique. 

Ihe actual value of vacuum thrust is given by equation 92 . The 
only measurement available for the SRB is the total pressure at the 
forward head end of the motor case and is modeled as 
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where b^ ^ and s^ ^ represents a bias and measurement noise respectively. 

The independent parameters, Aa and Ac^, are included in the 
model as states. The model assumed can be as given by equation 100 or 
another suitable dynamical process, i.e., random constant. 

The thrust deviation is given by the following truncated Taylor 
series as a function of the candidate independent parameters. 
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The partial derivatives for the v_ and m equations with respect to 
the independent pairameters are obtained directly from equation 106 . The 
mass deviation equation for the SRB is given as 
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The head pressure deviation > AP , is expanded sJjnilarly 
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3.3.3 Vehicle Mass 


The total rate of change of vehicle mass is given by 


T^(m ) — m__.,_ 
dt SSME 


"'SRB ^SSME ^SRB ''' "‘nON-CONSUME 

nom nom 


The first two terras in this equation are the apriori assumed nominal 
values. The third and fourth terms were discussed earlier. Ihe last 
term should be zero; however j it can include a mass bias uncertainty Amj^ 
Ihie equations, state and measurement, in which mass occurs are the 

•(B) 

V and equations. Assuming equation 110 C5n be summarized as 

• • 

m +• Amj^ then the mass can be written as ra + • Replacing this 

expression for the mass in the two indicated equations yields the 
following partial derivatives with respect to the Amj^. 
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4.0 PROJECTED ACTIVITIES DURING UPCOMING MONTH 


EXiring the next period, the filter portion of the estimation algorithm 
should bo checked out. Additionally, the smoother portion routines should 
have been coded. The smoother routines are more straightforward and represent 
less developments than the filter routines. 

"Rxe ability to bring up the NASA SRB predictive model on our computer 
will be determiend. If not successful, two other possibilities exist. 

The first, and least desirable, is to utilize the computer facility at the 
University of Florida under NASA sponsorship. Tlie second is to have pro- 
vided the necessary partial deriv^/.tives or allow SDI personnel access to 
these programs resident on the NASA computers for the purposes of generating 
the partial derivatives. There may be other possibilities than those 


mentioned here. 
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APPENDIX A 


PARTIAL DERIVATIVE OP THE VECTOR v Wt 0 


This partial derivative is one of several that occurs frequently in 
the formulation of the linearized system state and measuremotn egautions. 
The desired partial derivative is 


(cosecoet)v, + (Binvsinecost cos(pBin9coa+ ^ 

1 ~cos(psim|) 2 +sin(psxn(|> 


L. 

30 


(cos6sin4.)v, + (^InfslnOsln^, oos»einesin+, 

1 +cos<pcos(|» Vg -sxnqjcosfj' 3 


(-sin0)v^ H' (sin<pcos0 )v 2 + (cosfcos0) v^ 


A-1 


The resulting matrix is given in Table A-1 . 


TABLE A-1 . Partial Derivative of '*'C v wrt 


CD 


I 






cv 

fv 



> 



> 





A 

>*>> 


■a* 

■a* 


-a- 

•a- 


w 

1/5 


C4 

c< 


0 

P 


•H 

•H 


0 

0 


10 

(0 


©• 

e* 


e> 

£>• 


m 

H 


10 

P 


o 

*H 


Q 

*H 


u 

W 


0 

W 


+ 

1 


1 

+ 


■a- 



•a* 

a* 


c 

a 


w 

(0 



'H 

r- 

0 

0 

> 

W 

m 

> 

u 

u 


® 

<D 


® 

ICD 

•a- 

S 

g 

•a- 

c 

g 

c 

■w 

‘W 

(0 

>H 

'H 

■H 

w 

U1 

0 

(0 

(0 

W 

©■ 

e- 

u 

B- 

&• 

® 


u> 

ffl 

C 

IQ 

U] 

*rl 

0 

U] 

<H 

O 

<? 

V) 

u 

0 

(0 

0 

V 



o 



1 

1 

1 


+ 

+ 


ORIGIWAL PAG2 
OF POOR QUALO’Y 


•a* 

w 

o 

0 

CD 

C 

•H 

W 


CM 



CM 

(T 

> 

> 


> 

> 

a* 

•a- 


•a- 

■a- 

(0 

(0 


c 

0 

0 

. 

> 

*rl 


u 

0 


CO 

to 

® 

® 

•a- 

® 

CD 

(0 

to 

,c 

to 

to 

0 

0 

•H 

0 

o 

0 

u 

to 

u 

u 

e» 

e* 

m 

s- 

a- 

P 

to 

c 

p 

to 

•H 

o 

•H 

•H 

o 

¥ 


to 

1 




> > 



ffi 

<D 

r* 

g 

P 

> 


'H 


to 

to 

CD 

& 

©• 

CO 

c 

to 

O 

<H 

0 

O 

to 

0 




. 00 


■a" 

c 

•H 

CQ 

e- 

a 

•ri 

W 


•a- 

w 

o 

u 

CD 

.s 

m 

a- 

w 

0 

u 


■a- 

•S 

m 

e* 

CQ 

O 

u 


•a- 

03 

o 

o 

CD 

m 

e- 

c 

•H 

to 


CM 

no 



> 

> 







•a* 

•a* 



to 

to 



0 

0 



CJ 

u 



0- 

o 



g 

to 



•H 

0 



to 

0 



1 

+ 



-a- 

•a- 




g 

CM 

ro 

•H 

to 

•H 

to 

> 

> 

<D 

® 

CD 

CD 


c 

CO 

to 

•H 

•rt 

0 

0 

(0 

to 

0 

0 

€>• 

©• 

e- 

e- 

CO 

g 

10 

c 

0 

‘H 

0 

*H 

0 

to 

0 

to 

1 


ORIGINAL PAG^ fSf 
OP POOR QUALITY 


APPE^JDIX B 


av. 


m 


Ov 


0a 


(B)' 


15 V 


(B)' 


0v 


8B 

(B)' 


Qg 

0h' 


0S 

0h 


9v. 


and 


m 

0h 


Express ions 
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The resulting matrix is given in Table C-1 . 
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?iPPENDIX D 


AERODYNAMIC MODELING REGRESSION ANALYSIS AND RESULTS 


The aerodynamic data tables provided as IVBC3 data has been incorporated 
into an aerodynamic coefficient polynomial model. This modeling effort 
reduces the dimensionality of the numerical tables to one and reduces the 
storage requirements for the aerodynamic model. 

The coefficient model used for the two stages differ slightly as a 
result of the available data. The regression analysis led in the selection 
of the form of the aerodynamic model. Terms with insignificant correlation 
were eliminated from the model . 

In equation form, the first stage static coefficients of axial force, 

C^; normal force, Cj^; pitching moment, C^; rolling moment, C^} side force, 

C^; and yawing moment, C^; are given below 
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The corresponding second stage model is given by the following equations 
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For the first stage, data from an angle -of -attack range of -6 to +6 degrees 
was used in the regression analysis. Data from a range of -8 to +4 degrees 
was used for the second stage. The results, TCXX..., from the regression 
analysis is presented below for each of the coefficients, C 


, above . 
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table D— 2 . First Stage Normal Force Coefficient 
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